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The free energy between a static quark and an antiquark is studied by using the color-singlet 
Polyakov-line correlation at finite temperature. We perform simulations on 32 3 x 12, 10, 8, 6, 4 
lattices in the high temperature phase with the RG-improved gluon action and 2+1 flavors of the 
clover-improved Wilson quark action. Since the simulations are based on the fixed scale approach 
that the temperature can be varied without changing the spatial volume and renormalization factor, 
it is possible to investigate temperature dependence of the heavy-quark free energy without any 
adjustment of the overall constant. We find that, the heavy-quark free energies at short distance 
converge to the heavy-quark potential evaluated from the Wilson-loop operator at zero tempera- 
ture, in accordance with the expected insensitivity of short distance physics to the temperature. 
At long distance, the heavy-quark free energies approach to twice the single-quark free energies, 
implying that the interaction between heavy quarks is screened. The Debye screening mass ob- 
tained from the long range behavior of the heavy-quark free energy is compared with results of 
the thermal perturbation theory and those of Nf = 2 and Nf = lattice simulations. 
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1. Introduction 

An interaction between a heavy quark and an antiquark in hot QCD medium is one of the 
most important quantities to characterize internal properties of the quark-gluon plasma (QGP). 
Experimentally, it is related to the fate of the charmoniums and bottomoniums in QGP created in 
relativistic heavy ion collisions. The interaction between heavy quarks at finite temperature (T) 
have been studied via heavy-quark free energies in the quenched approximation [2|, [3]], in Nf = 2 
QCD with the staggered [|] and the Wilson quark actions [g, g], and in Nf = 2 + 1 QCD with 
the staggered quark action ^ In this report, we present resent studies on the heavy-quark free 
energy in dynamical simulations of Nf = 2 + 1 QCD with the Wilson quark action. We perform 
finite-temperature simulations ^ on 32 3 x (12-4) lattices adopting the fixed scale approach which 
enables direct comparison of the heavy-quark free energy at different T since the spatial volume 



and the renormalization factor are common to any T [JlOp . The corresponding zero-temperature 
configurations are taken from the results of the CP-PACS and JLQCD Collaborations [O]. The 
free energies between heavy quarks at finite T are evaluated from the Polyakov-line correlator and 
compared with a heavy-quark potential at T = calculated from the Wilson loop operator. We also 
extract the Debye screening mass, mo{T), by fitting the heavy-quark free energy with the screened 
Coulomb form, and compare it with results of thermal perturbation theory and of lattice simulations 
in Nf = 2 and N f = QCD. 

It is found that the heavy-quark free energy is insensitive to T at short distance region, whereas 
it is screened by thermal medium at long distance region and converges to twice a single-quark free 
energy. The Debye screening mass shows better agreement with that calculated in the next-to- 
leading order of thermal perturbation than in the leading order. Also it receives a sizable effects 
from the dynamical light quarks. 



2. Heavy-quark free energy 

At T = 0, interaction between a static quark and an antiquark can be studied by the heavy- 
quark potential evaluated from the Wilson loop operator. The resulting potential takes the Coulomb 
form at short distances due to perturbative gluon exchange, while it takes the linear form at long 
distances due to confinement: 

V(r) = - — + ar+V . (2.1) 

For T > 0, inter-quark interaction may be studied by the heavy-quark free energy F 1 (r,T) evaluated 
from a Polyakov-line correlation function in the color-singlet channel with Coulomb gauge fixing 
[0: 

F\r,T) = -rin(Trn t (x)n(y)), (2.2) 

where r = |x — y| and Q(x) = Yl^' =l t/4(x, x) with U^(x, x) being the link variable in the temporal 
direction. At zero temperature, we expect F 1 (r,T = 0) = V(r), while at high temperature, F 1 (r,T) 
is screened by thermal medium and may behave as the screened Coulomb form: 

F \rj) = -^H e -™o(T)r ^ (23) 
r 

where a and mrj are the effective coupling and the Debye screening mass, respectively. 
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3. Fixed scale approach 

In the conventional fixed N t approach where T is varied by changing the lattice spacing a, V(r) 
receive different renormalization at each T, and thus are usually adjusted by hand such 
that V(r) and F (r, T) coincide with each other at a short distance, assuming that the short distance 
properties are insensitive to the temperature. On the other hand, in the fixed scale approach, the 
temperature T = (N t a)~ l is varied by changing the temporal lattice size N t at fixed a [|T0|]. 

In the fixed scale approach, because the coupling parameters are the same for all temperatures, 
the renormalization factors are common to all temperatures. The spatial volume is also the same. 
We can thus directly compare the free energies at different temperatures without any adjustment. 
We show below that F l (r,T) for different T approaches to V(r) at short distances, which proves 
the expected insensitivity of the short distance physics to the temperature. 



4. Results of the lattice simulations 

We employ the renormalization-group improved gluon action and 2 + 1 flavors of nonpertur- 
batively 0(a) -improved Wilson quark actions. Zero-temperature configurations are given by the 



CP-PACS and JLQCD Collaborations [JTT]] at m n /m p = 0.6337(38) and m K /m K * = 0.7377(28). 
Finite temperature simulations with the same parameters are performed on 32 3 x N t lattices with 
N, = 12, 10, 8, 6 and 4, which correspond to T ~ 200-700 MeV |g|. We generate full QCD con- 
figurations by the hybrid Monte Carlo algorithm and measure the heavy-quark free energy using 
500-700 configurations at every 5 trajectories after thermalization. The statistical errors are deter- 
mined by the jackknife method with the bin-size of 20 trajectories. The absolute scale is estimated 
from the Sommer parameter, ro = 0.5 fm. 

Figure shows the results of the heavy-quark potential V (r) at T = and the heavy-quark free 
energies F l (r,T) at various temperatures as functions of r. At T = 0, V(r) shows Coulomb-like 



and linear-like behaviors at short and long distances, respectively. A fit with Eq. (2.1 ), as shown by 
the dashed line in Fig. [I], gives a = 0.441 and y/a = 0.434 GeV. 

For T > 0, we note that the heavy-quark free energies F x (r, T) at all temperatures converge to 
V(r) at short distances. This means that the short distance physics is insensitive to temperature. As 
stressed above, unlike the case of the conventional fixed-A 7 , approach in which this insensitivity is 
assumed and used to adjust the constant terms of F l (r, T), our fixed scale approach enabled us to 
directly confirm this theoretical expectation. 

At large r, F l (r, T) departs from V(r) and eventually becomes flat due to Debye screening. In 
Fig. H, the asymptotic values of F 1 (r, J) at long distance are also compared with 2Fq denoted by 
the arrows where the thermal average of a single Polyakov-line is defined as Fq = — rin(Tr£2). We 
find that F l (r, T) converges to 2Fq quite accurately at long distances. 



5. Debye screening mass 

In order to extract the screening mass mo, we fit F l (r,T) — 2Fq by the screened Coulomb 
form, Eq. (2.3). To determine the appropriate fit range, we estimate the effective Debye mass from 
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Figure 1: Heavy-quark free energies at various temperatures. The heavy-quark potential at T = was 
calculated by the CP-PACS and JLQCD Collaborations [juj]. The fit result of V(r) by the Coulomb + linear 
form is shown by the dashed line. The arrows on the right side denote twice the thermal average of the 
single-quark free energy defined as 2Fq = — 27Tn(Tri2). 
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(5.1) 



' F l (r + Ar)-2F Q 

Investigating the plateau of nf^{T;r), we choose the fit range to be 0.35 fm < r < 0.56 fm. Results 
of mo{T)/T are shown in the left panel of Fig. ||, and does not have strong dependence on T. This 
property is qualitatively consistent with the prediction in thermal perturbation theory: mo ~ gT 
where g is the running coupling. 

To make a quantitative comparison to the result of the thermal perturbation theory, we define 
the 2-loop running coupling by 



S 21 2 (AO = /3oln(-^) 2 + f In 

A MS P0 



ln( 



A MS 



(5.2) 



with the QCD scale parameter A^T = 260 MeV Q13|], where we assume a degenerated Nf = 3 
case. The renormalization point }X is assumed to be in a range /J. = nT-3nT. Then, the Debye 
mass as a function T in the leading-order (LO) thermal perturbation theory is given by: 



(5.3) 
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Figure 2: Left: Results of the Debye mass for Nf = 2+1 obtained from the heavy-quark free energy with 
those of the leading order (LO) and next-to-leading order (NLO) thermal perturbation theory for Nf = 3. 
Right: Comparison among the Debye masses for Nf = 2 + 1, Nf = 2 [§] and N t = [[ToJ] extracted from the 
heavy-quark free energy on the lattice. Arrows in the right panel indicate critical temperatures in Nf = 2 and 
Nf = cases. 



where we have neglected the quark mass effects. Dashed lines in Fig. || (left) represent /n^ (r) for 
Nf = 3 for the above range of the renormalization point. We find that the LO Debye mass does not 
reproduce the lattice data. 

A formula in the next-to-leading-order (NLO) is also available from the resummed hard ther- 



mal loop calculation [14]: 



in 



NLO 



l + -Z-87l(T) 
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l mag 
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(5.4) 



where m mag (T) = C m g 2 (T)T is the magnetic screening mass. Since the factor C m cannot be de- 
termined in perturbation theory due to the infrared problem, we adopt C m ~ 0.482 calculated in a 
quenched lattice simulation [15] as a typical value. Solid lines in Fig. || (left) represent the NLO 
results for Nf = 3. We find that the NLO Debye mass is approximately 50 % larger than the LO 
Debye mass and shows a better agreement with the lattice data. 

Finally we study the flavor dependence of mo. Fig. ^ (right) shows trip in lattice simulations 
of Nf = 2 + 1 QCD (this work), that of Nf = 2 QCD with an improved Wilson-quark action at 
m n /m p = 0.65 and that of the quenched QCD (Nf = 0) JT0|]. Arrows on the horizontal axis 



indicate critical temperatures for Nf = 2 (T c ~ 186 MeV at m^/rtip = 0.65 [16]) and Nf = (T c ~ 
290 MeV). We find that nip for Nf = 2 + 1 is comparable to that for Nf = 2, whereas it is larger 
than that for Nf = 0. A similar result was obtained with a staggered-type quark action [j|]. 



6. Summary 

We studied the free energy between a static quark and an antiquark at finite temperature in 
lattice QCD with 2+1 flavors of improved Wilson quarks on 32 3 x (12-4) lattices in the high 
temperature phase. We adopted the fixed scale approach which enables us to directly compare 
the free-energies at different temperatures without any adjustment of the overall constant. At short 
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distances, the heavy-quark free energies, evaluated from the Polyakov-line correlations in the color- 
singlet channel, show universal Coulomb-like behavior common to that of the heavy-quark poten- 
tial at zero temperature evaluated from the Wilson-loop operator. This is in accordance with the 
expected insensitivity of short distance physics to the temperature. At long distances, the heavy- 
quark free energies approach to twice the single-quark free energies calculated from the thermal 
average of a Polyakov loop. Also, we extracted the Debye screening mass mj)(T) from the heavy- 
quark free energy and found that the next-to-leading order thermal perturbation theory is required 
to explain the magnitude of mj)(T) on the lattice. Comparison to the previous results at Nf = 2 and 
Nf = 0, shows that the dynamical light quarks have sizable effects on the value of mj){T). 
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